Dynamics of Fractures in Quenched Disordered Media. 
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We introduce a model for fractures in quenched disordered media. This model has a deterministic 
extremal dynamics, driven by the energy function of a network of springs (Born Hamiltonian) . The 
breakdown is the result of the cooperation between the external field and the quenched disorder. 
This model can be considered as describing the low temperature limit for crack propagation in solids. 
To describe the memory effects in this dynamics, and then to study the resistance properties of 
the system we realized some numerical simulations of the model. The model exhibits interesting 
geometric and dynamical properties, with a strong reduction of the fractal dimension of the clusters 
and of their backbone, with respect to the case in which thermal fluctuations dominate. This 
result can be explained by a recently introduced theoretical tool as a screening enhancement due 
to memory effects induced by the quenched disorder. 
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I. INTRODUCTION 

In the last years, many models have been proposed 
to describe the formation of cracks in different kind of 
materials They are inspired by the study of non 

equilibrium fractal growth processes, like the Dielectric 
Breakdown Model (DBM) || and Diffusion Limited Ag- 
gregation (DLA) |(| . These models are based on two dif- 
ferent mechanisms for the fracture growth: 1) A stochas- 
tic effect due to the thermal fluctuations in the medium, 
driven by an external field. 2) A deterministic dynam- 
ics, when the main source of randomness is the quenched 
disorder of the medium (like in the Invasion Percolation 
(IP) model 0. ), being the thermal fluctuations negligi- 
ble (i.e. low temperature limit). 

A very interesting model, belonging to the first class, 
is the Born Model (BM) [|]|| , where minimization of the 
elastic energy is used to compute the driving field. This 
model for fracture propagation is the analogous of the 
DBM for Laplacian growth. It describes the fractures 
propagation as a stochastic process, where the proba- 
bilistic mechanism represents growth instabilities, like for 
example density fluctuations in a gas. In this approach 
the quenched disorder of the medium is negligible. 

In this paper we want to consider the limit of low tem- 
perature for the BM, where the driving field cooperates 
with quenched disorder to produce the breakdown pat- 
terns. The quenched disorder can be thought to represent 
the effect, at a mesoscopic scale, of defects of the break- 
ing layer. In this version of the model, at each time step 
the growing bond is selected deterministically, i.e. the 
bond with the smallest ratio between the local driving 
field and the quenched disorder grows. The system is a 
2d triangular lattice of springs. We apply to two parallel 
boundaries of the system an uniaxial tension and fixed 
boundary condition to the others. In this way we ob- 
tain directed, crack patterns, orthogonal to the applied 



stress. The system has two independent length scales, 
the height h and the width L. This allow us to evaluate 
more clearly the various scaling regimes of these struc- 
tures. Usually, with this boundary conditions, two differ- 
ent phases, called respectively scaling regime and steady 
state, are present |fO| . 

We will see that in quenched version of BM (QBM), be- 
cause of the very strong screening on the growth process, 
the first regime (where many different branches compete 
during the growth) is nearly absent. In fact, most of the 
growth dynamics develops during the steady state, where 
only one branch survives and the cluster is statistically 
self-similar with well defined fractal properties. Together 
with the fractal dimension, the backbone and chemical 
distance exponents characterize completely the fractures 
produced by our model. 

We compare our results with those obtained previously 
for the stochastic version of the BM || showing that the 
fractal dimension of the clusters and of the backbone, in 
similar conditions, are strongly reduced. We explain this 
as a consequence of the absence of thermal fluctuations; 
the quenched disorder produces memory effects giving 
screening effects similar to those of IP [jl2| . On the other 
hand, in the same conditions, the screening effect pro- 
duced by the modulation of the driving field gives a frac- 
tal dimensions much smaller than IP, which is the limit 
of the model in which there is no external driving field 
or it is constant. The same qualitative result has been 
recently found, both numerically and analytically, for the 
quenched version of the DBM model [Q. Therefore, the 
numerical results found here for the QBM support the 
feeling that this screening enhancement is a general prop- 
erty of all deterministic models with quenched disorder 
and a driving field. 

The distribution of the quenched disorder related to 
grown (broken) sites (acceptance profile) is also studied 
for different values of the model parameters. This dis- 
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tribution reaches, during the evolution, an asymptotic 
shape. In the case of Invasion Percolation, the asymp- 
totic distribution is a step function, with the discontinu- 
ity at the value of quenched disorder coincident with the 
critical point of classic percolation QJl^l- This step-like 
shape of the acceptation profile indicates that the dynam- 
ics develops avalanches with a scale invariant distribution 

In our model, the presence of the stress field, how- 
ever, does not allow the presence of a critical threshold 
in this asymptotic distribution of quenched disorder. As 
a consequence, the dynamics does not develop scale in- 
variant avalanches, but these avalanches have a typical 
size fl2] , |l3|| . A very important universality relation can 
be written, which explains, in terms of the dependance of 
the fractal properties on the parameters of the model, the 
cooperation between the driving field and the quenched 
disorder in developing such a fractal structure. 

The paper is organized as follows. In section II the 
model is introduced and compared with the correspond- 
ing stochastic model. The details for the realization of 
simulations are specified. In section III we describe our 
numerical results for the fractal properties of the clus- 
ters, the backbone and the chemical distance, for the 
roughness exponent of the chemical distance, and for the 
statistical effective distribution of quenched variables on 
the growth interface. In section IV the universality rela- 
tion is demonstrated and a theoretical tool is introduced 
to explain the two screening effects in producing fractal 
properties. Finally in section V we discuss the results 
and draw some conclusions. 

II. THE MODEL 

The BM describes the medium as a discrete set of 
springs. The equilibrium state is obtained by imposing 
minimization of the energy while dynamics of fracture 
is given by assigning a rule of growth. For what con- 
cerns the equilibrium state, we imagine to break only 
one spring a time to model a system of slowly developing 
fractures, and after every breakdown a new equilibrium 
state is computed. We describe the energy of the sys- 
tem by means of the same potential energy. This elastic 
potential energy consists of two different terms, describ- 
ing, respectively, a central force and a non-central force 
contribution: 

= 2 Y,( a - ^K* - u o) ■ nj} 2 + P[ui - uj} 2 , (i) 

where Ui is the displacement vector for site i, fij is the 
unity vector between i and j, a and (3 are force constants, 
and the sum is over the nearest neighbor sites connected 
by a non-broken spring. 



For any equilibrium state must result V^} V(ui, Uj) = 
0,Vi, where j n.n. of i. This condition yields to a series 
of equation for the Hi's that can be solved imposing the 
boundary condition. The initial boundary condition is a 
uniform dilation on the left and right side of the sample. 
This boundary condition change, taking into account all 
the springs broken during the evolution of the crack. 

For what concerns the rule of growth we have chosen 
a deterministic rule selecting the bonds to break accord- 

1/2 

ing to its "generalized elongation" V i j . The particular 
rule explained in the following equations has been in- 
spired by the analogies with the DBM model. We think 

1/2 

V i j as a field acting on the spring between sites i and 
j. Since the system is characterized by the presence of 
random quenched defects (represented as a quenched ran- 
dom noise), we assign to each spring a random number 
Xij extracted by the probability density: 

Po(sij) = Z^-^J 1 . ( 2 ) 

^ sup 

where the parameter a (a £ [0, oo)) modulates the im- 
portance of the disorder in the mechanical properties of 
the system, and the variables are defined in the range 
[0,x sup ]. From eq. ^ one can derive the mean value (x) 
of the disorder assigned to the bonds: 

(x) = / dxxp (x) = — — x sup . (3) 
Jo a + 1 

A "fragile" material correspond to small values of a 
((x) ~ 0), while a "rigid" one corresponds to big val- 
ues of a ((x) ~ x sup ). Then we define a set of dynamical 
variables: 

V l ,j(t) = A hJ (t)x hJ (4) 

where Aij(t) = y -(t) 1 / 2 • Then, at each time step, the 
spring with the smallest value yij(t) on the growth in- 
terface of the fracture is broken. 

It is customary to modulate the influence of the stress 
field on the growth dynamics with the introduction of a 
parameter 77. In this case the formula for Aij becomes 
= v (t)'?/ 2 • When rj — there is no field and 
the model has the same dynamics as invasion percolation 
, while rj — 00 correspond to an infinite strength of 
the field and the cluster is an one-dimensional straight 
line. In fact, the disorder effects are negligible when 
compared to an infinite field. Some details have to be 
specified since in principle their variation could affect the 
fractal properties H^PI. Due to the vectorial nature 
of eq.(Q) a triangular lattice is more appropriate to model 
the medium. In fact, for a squared lattice when j3 = 0, 
the system behaves as a set of independent planes with- 
out connection between each other. For this reason we 
will follow ^,0, by considering only triangular lattices. 
Furthermore, the growth interface at any time is given 
by the set of non broken bonds nearest neighbors to the 
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cluster ones (i.e. the perimeter of the fracture cluster). 
This corresponds to the implementation of a connectivity 
condition (FIG. 0). 

III. NUMERICAL SIMULATIONS 

We performed several realizations on systems of size 
LxL (L = 64, 128), in triangular geometry (with periodic 
boundary conditions on top and bottom sides). 

A quenched variable Xj t j is assigned to each bond 
The Xij 's follow the distribution function of eq. (|^) with 
x S up — 0.5. In fact, a constant strain equal to 0.5 lattice 
units is applied in the horizontal direction (see FIGj|). 
At each time step, the stress field over the interface bonds 
is computed. Then the bond with the smallest value 
Di.j is broken, the new stress fields are computed and 
new bonds are added to the perimeter. By iterating this 
dynamics, one obtains structures like that shown in FIG. 

The most important quantity characterizing the struc- 
tures generated by the model is the fractal dimension 
Df. For each value of L, we have generated a set of 
20 realizations and computed the fractal dimension by 
the box-counting method. This has been done for differ- 
ent values of P/a. In fact, in the equilibrium condition 
V {g.yV (Hi) = we deal only with the ratio (3/ a. For this 
reason, we decided to vary f3 and to keep a — 1. Further- 
more, we performed the same analysis also by varying r\ 
(the parameter modulating the effects of the field) and 
a (the parameter modulating the disorder). Our results 
are shown in Tab. § and Tab. || [n]. 

The model shows a continuous dependence of the frac- 
tal dimension on the parameter (3, as found in || for the 
stochastic BM. The dependance of the fractal dimension 
on the parameters rj and a is, instead, new and interest- 
ing. In fact, from our simulations we see that the fractal 
dimension actually depends only on the product an. This 
introduces a precise relationship between the indices de- 
scribing the properties of the medium and the properties 
of the stress field. 

If we compare this result with the case of invasion per- 
colation, where the field is absent and the fractal dimen- 
sion does not depend on the value of a Jj],[l2| , we see that 
the introduction of the stress field breaks the symmetry 
with respect to a, leading to less universal "critical" prop- 
erties. In the following section we present an analytical 
demonstration of this universal property. 

As it can be noted, for fixed a, the fractal dimension 
of the fracture cluster, for any value of n > 0, is less 
than the IP one. This is due to the fact that for n > 0, 
there is a screening effect related to the physical stress 
field, in addition to the screening related to the mem- 
ory effects of the quenched disorder. Moreover, our nu- 
merical results show clearly that in a "fragile" material 
(small a, (x) ~ 0) the fractures have a big fractal di- 
mension, while in a "rigid" material (big a, (x) ~ x sup ) 



the fractures tend to be straight lines with fractal di- 
mension near to 1. This result sounds qualitatively rea- 
sonable from an experimental point of view. A "fragile" 
material could correspond to a material with many im- 
purities which lower its resistance to rupture, allowing 
many bonds to be broken. A rigid material could repre- 
sent a material without impurities at zero temperature, 
in which fractures are straight lines along the direction 
of maximum stress. However, a quantitative comparison 
with experiments is still not accessible, since it needs a 
clear connection between what we call quenched disorder 
in our model, which we believe to give a description of 
the system at a mesoscopic scale, and the microscopic 
disorder in real systems. 

In the next section we will propose an analytical ex- 
planation for our numerical findings. The same qualita- 
tive results have been recently found for a similar model, 
the quenched dielectric breakdown model (QDBM) ]l3) . 
This suggests that the dependance on the product of a 
and r\ of the dynamics and geometry of the patterns plus 
a screening enhancement effect are general properties of 
all the deterministic models with driving field in presence 
of quenched disorder. 

A further characterization of the topological and con- 
nectivity properties of the aggregates is given by the ex- 
ponents ruling the scaling of two subsets of the clusters: 
the chemical distance and the backbone. The chemical 
distance is the shortest path between the two ends of 
the aggregate and shows interesting self-affine properties. 
The backbone is obtained cutting from the cluster all the 
tips and the dangling loops connected to the chemical dis- 
tance (see FIG. ||). This part of the cluster influences the 
macroscopic transport properties of the medium, while 
the chemical distance gives the shape of the line separat- 
ing the system in two parts after the breakdown. It is 
worth to point out that in models with a scalar external 
field, like DBM §] and QDBM H), no closed loops are 
allowed, and the backbone trivially coincides with the 
chemical distance ||. 

To determine the backbone and the chemical distance 
we followed a method used in || , based on the topological 
properties of the clusters. As an example, in FIGG. |^, 
U we show, respectively, the backbone and the chemical 
distance we obtain for the cluster of FIG. ||. In Tab. |^ 
we collect the results of a box-counting analysis. The 
backbone fractal dimension is significantly larger than 
one, although quite smaller than for the stochastic Born 
model M. Instead, the fractal dimension of the chemical 
distance is 1 for both L = 64 and L = 128. 

We then studied the roughness exponent x of the chem- 
ical distance, which gives the scaling of the mean squared 
lateral width W(l) of a self-affine path with respect to its 
length I. In our case, the path develops vertically, and its 
length is measured along the vertical direction, while the 
width is computed along the horizontal direction. Then, 
if is the ^-coordinate of a point on a chemical dis- 
tance whose y-coordinate is j, W(l) is defined as: 
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W(l) 



jo+l 
3=3o 



(5) 



where (...) represents a mean over all portions of length I 
of the chemical distance and over different realizations of 
quenched disorder, and i is the mean horizontal position. 
For a self-affine path W(l) ~ l x holds. In FIG. @ we 
show the scaling behaviour of W 2 (l) versus I, for clusters 
of size L — 128 and for different values of the parameter 
rj. A least square fit, reported in the figure, gives for 
values of rj ranging over an order of magnitude (from 0.2 
to 3), values which seem quite stable within error bars. 
For rj = 3 we get x — 0-62 ± 0.03, for rj = 1.0 we get 
X = 0.64 ± 0.03 and for rj = 0.2 we get x = 0.66 ± 0.03. 

The independence of x from rj (i.e. from the product 
arj) indicates a very important universality property of 
the model. In particular the values we find are quite sim- 
ilar to that of the directed percolation (x — — 0.625... 

[p~9|| ) , so we can suppose that, for what concerns the 
chemical distance, the model is in the same universality 
class of directed percolation. 

The evolution of the dynamics towards an asymp- 
totic stationary state is characterized by the acceptance 
profile a(x). It gives the rate of acceptation (selection 
by the growth process) for the values of the quenched 
variables associated to the grown bonds, in the range 
[x,x + dx]. We have performed a set of 10 realizations 
of size 512 x 512, for r\ = 1.0, 2.0, 3.0 and a = 1, \3 = 0.5. 
This allowed us to follow the time evolution of a(x) up 
to t = 500, where the asymptotic state is reached. In 
FIG. H we show the final shape of a(x), for the different 
values of rj. The presence of the stress field prevents the 
dynamics from developing scale invariant avalanches and 
this is reflected in the absence of a critical threshold in 
a(x) [fl2"|| . However, as the value of rj becomes smaller, 
and the role of the stress field becomes less relevant, one 
can see that a(x) seems to develop a discontinuity. What 
we expect is that in the limit rj — one recovers the IP 
dynamics. Also in this case we obtain results that agree 
with those obtained for the QDBM model |l]|. In fact, a 
common picture, that we describe in a longer paper, can 
be used to study both models [BOfl. 



IV. THEORETICAL RESULTS 

In this section we present some theoretical results con- 
cerning the universality properties of the model and the 
explanation of the screening effects which give rise to such 
fractal breakdown patterns. 

To demonstrate that all the dynamical and geometrical 
properties of the model depends only on the product arj 
instead on the two parameters separately, consider the 
generic variable 



Vi,j (t) 



(6) 



where extracted from the density function 

^ sup 



(7) 



Now perform the following transformation of the 
quenched variables 

*ij = KA a (8) 

Introducing eq.(||) in eq.(0) one can note that the 
variables Eire uniformly distributed between and 
z S up = Xsup- So the density function of the new vari- 
ables does not depend on the parameter a apart from 
the value of z sup (unique for all the bonds) value. Intro- 
ducing eq.(||) in eq.(||) we obtain 



K,,i(t)(»* 2 



l/a 



(9) 



Finally, as the dynamics is extremal it does not change if 
instead of the variables j/j t j (t) we consider as bond vari- 
ables 



(10) 



This equation shows that all the dynamical and ge- 
ometrical properties of the model depends only on the 
product arj as asserted before. 

Now we switch to the explanation of how the two 
screening effects (the geometrical one and the proba- 
bilistic one) cooperate in the model during the dynam- 
ical evolution. In doing so we will briefly introduce a 
generalized version of the Run Time Statistics (RTS), 
which was first introduced to study extremal dynamics 
in quenched disordered media, like IP ( E3, p2|). The 
RTS approach consists in a transformation of the deter- 
ministic quenched extremal dynamics in a stochastic one 
through the introduction of effective, time-dependent, 
density functions for the dynamical variables y^j. The 
role of these densities is to store information on the past 
growth history of the system (ljjl^l . 

This approach allows us to have a well defined growth 
probability distribution at any time step and an updating 
rule of the effective density functions after any elemen- 
tary growth event. 

To explain simply the observed screening effects, sup- 
pose to have only two bonds, 1 and 2, in the interface 
at time to an d consider a — 1 and r) > 0. Let the 
quenched variables of the two bonds be respectively X\ 
and X2, which are uniformly distributed between and 
1. Moreover let the related stress fields be respectively 
Ei j = V"i(t) r '/ 2 and E 2 , to = V 2 {t)' n/ ' 2 where for example, 
without loss of generality, E\ >to > Ei^ - Then we have 
Vi.ta = Xi/Ei t t for i = l,2 uniformly distributed between 
and l/E itto , i.e. 



Pi,t (y) = E itto i = l,2 



(11) 
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Pi,t (y) represents the effective density of the bond i at 
time to- We want now to calculate the probability /ii,t 
that bond 1 grows at to + 1, i.e. that 2/1, t < V2,t - This 
probability is given by: 



1/-Bi,t 



1/-B2,t 



d yi E Uo l d y 2 E 2M =l-^^>\ (12) 



Now, we suppose this is the real growth event and com- 
pute the new probability density p 2 ,t a +i{y2) of the vari- 
able j/2,t +i conditioned to that event. For sake of sim- 
plicity, we assume that the stress field of the bond 2 keeps 
constant, i.e. 2/2,t +i = 2/2, t - By applying the rules of 
conditional probability, we obtain: 

P2,t +i(y2)=P2,ro(y2) / dyiPi,t {y x )6(l/E 1M -yi) (13) 
Jo 

where the function 9 is the step function and indicates 
that the variable j/i must be smaller than 1/Eit . By 
introducing eq. O in eq. [13 we obtain: 



P2,t +i(V2) = { E 2 T' t 



%^y2;0< y2 <^ 



Mi, t ' -Ei, to 



1 <2/2< 1 



(14) 



E2. t 



This function is represented in FIG. ||-b and, the effective 
probability density of yi at time to + 1 is more concen- 
trated towards the high values of the variable than at 
time to at which was uniform. If we introduce this den- 
sity in the eq. [l^, we see that the growth probability 
of bond 1 increases, and by normalization the growth 
probability of bond 2 decreases. This is what is called 
before " probabilistic screening effect" , which is nothing 
else than a memory effect. In deriving the eq.([l4|), we 
made the hypothesis that the stress field of the bond 
2 remains constant. In general this is not the case, in 
fact when a bond i grows at time t the stress field of a 
bond j remaining on the growth-interface decreases, i.e. 
Ej } t+i < Ej >t . This leads to the second kind of screening: 
the geometric one. In fact, supposing E^to+i < -^2,t , we 
obtain: 



2/2,i + l 



E 2,t 



E- 



2,to + l 



-£/2,t ^ ^2,to + l2/2,t 



(15) 



Note that, as ^2,^+1 < 1, we have y2,t +i < 2/2, t - From 
cqs.(|i"4|) and ( [i"5| ) we obtain the effective probability den- 
sity of the dynamical variable of the bond 2: 



P2,t +l(V2) 



. r Mi, t ^2,t + i ' — »2 _ El tg 
' X E 2 ,t + i . ^2,tp + l < < 1 

Mi,t ' -Ei, t0 - y2 - E2,t + 1 



-42,tp + l 



(16) 



where now we account for the decrease of the stress field 
in the variable yi- Since, by definition, Ai^+i > 1, this 
function is more concentrated towards the high values of 
the variable than in the previous case (Fig. ^-c). Thus, 
the growth probability of bond 2 further decreases with 



respect to eq.jl^). This reflects the presence of both the 
screening effects, the geometric one due to the decreasing 
stress fields of "old" perimeter bonds, and the temporal 
one produced by memory effects. 

In this way we have shown simply how the two kinds 
of screening effects add the evolution of the dynamics 
leading to fractal structure with small fractal dimensions. 

The generalization of these arguments to the case in 
which the interface is composed by many bonds with dif- 
ferent initial probability densities, leads to the formula- 
tion of the RTS for extremal dynamics with quenched 
disorder and a driving field, which has been introduced 
in jl3| to study the problem of the electric discharge in 
disordered dielectric system. 



V. CONCLUSIONS 

We have introduced and studied a model for the 
dynamics of fractures in a quenched random medium. 
The model exhibits interesting properties, for example a 
strong reduction of its fractal dimension with respect to 
the case of a dominating thermal noise (stochastic dy- 
namics H), and a roughness exponent for the chemical 
distance whose value seems to be independent on the 
model parameters. Moreover, the fractal dimension of 
the fractures depends continuously on the material prop- 
erties, represented in the model by a power law distri- 
bution of quenched disorder with a tunable parameter a. 
This is, from an experimental point of view, reasonable. 
In a situation in which thermal fluctuations of the lat- 
tice are irrelevant, the only source of noise, and thus of 
fractal properties, in the growth process is the intrinsic 
quenched disorder of the samples (impurities, vacancies 
etc.). Our model, however, is still too an idealized one to 
allow a quantitative comparison with experiments. 

It would be interesting to compute analytically the 
fractal dimension of this model by using the Fixed 
Scale Transformation (FST ^l|) approach, after having 
mapped its dynamics onto a stochastic one, with the RTS 
method. The mapping is in fact possible, but the approx- 
imations one has to use in the FST approach give rise to 
quite poor numerical results, as reported in a longer pa- 
per on this subject [^(J . 

The model can be further generalized. A research di- 
rection we are presently following concerns the unifica- 
tion of the low temperature (deterministic, extremal dy- 
namics with quenched disorder) and high temperature 
(stochastic dynamics with thermal noise) regimes, by the 
introduction of a temperature-like parameter T which 
can tunc the transition between the two regimes p2~| . 
This could represent a further step toward a more real- 
istic description of fracture propagation in solids. One 
interesting property of this generalized model is that in 
the high temperature limit, the level of approximation 
of the FST method becomes excellent pj[ . This could 
allow us to get a good estimation of the low temperature 
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fractal dimension by an extrapolation of the values we 
get in the high temperature case. 



TABLE I. Behavior of the fractal dimension of the QBM 
model for different values of /3, with a = 1,tj = 1, for sizes 
L = 64, 128. 
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p 


Df (L = 64) 


D f (L = 128) 


0.0 


1.15 ±0.03 


1.14 ±0.02 


0.005 


1.18 ±0.03 


1.16 ±0.02 


0.05 


1.20 ±0.03 


1.17 ±0.02 


0.5 


1.22 ±0.02 


1.20 ±0.02 


5 


1.26 ± 0.02 


1.24 ±0.02 



TABLE II. Fractal dimension of the QBM model for differ- 
ent values of rj and a for clusters of size tL = 64. The values 
of the other parameters are: /3 = 0.5 and a = 1. 



V 


£>/(!,, a =1) 


a 


D f (r) = l,a) 


0.5 


1.36 ±0.03 


0.5 


1.33 ±0.03 


1.0 


1.21 ±0.03 


1.0 


1.21 ±0.03 


2.0 


1.13 ±0.03 


2.0 


1.15 ±0.03 


3.0 


1.10 ±0.03 


3.0 


1.10 ±0.03 



TABLE III. Fractal dimension of the QBM model for dif- 
ferent values of n and a for clusters of size L = 128. The 
values of the other parameters are: (3 — 0.5 and a = 1 



1J 


D f {n,a=l) 


a 


£>/(»? = 1, a) 


0.5 


1.42 ±0.02 




0.5 


1.40 ±0.02 


1.0 


1.20 ±0.02 




1.0 


1.20 ±0.02 


2.0 


1.15 ±0.02 




2.0 


1.16 ±0.02 


3.0 


1.11 ±0.02 




3.0 


1.10 ±0.02 


TABLE IV. Behavior of the 


backbone fractal dimension 


{D B 


) and (Dc) of chemical distance fractal dimension, with 


P = 


0.5, a = 1, r) = 1, a = 


1, for clusters of size L — 64 and 


L = 


128. 








L 


D B 


D c 


64 


1.07 ±0.02 




1.00 ±0.02 


128 


1.10 ±0.02 




1.00 ±0.02 




FIG. 1. Connectivity condition. Broken bonds are indi- 
cated by dashed lines, while interface bonds are indicated by 
thick lines. 
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FIG. 2. Uniaxial tension applied to the lattice. 
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■ mean square width L-128, r|=3 
f(l)=a P, X=0.62, Az=0.03 

° mean square width L=128, r|=l 
g(])=b P, x=0.64, Ax=0.03 

» mean square width L-128, r|=0.2 

i^X n ti: i iiiii 



h(l)= 



", z=0.66, Ax=»-»3 



X. 



FIG. 3. A crack pattern generated by our model, with pa- 
rameters rj = 1.0, a — 1.0, a — 1.0, P — 0.5. 



FIG. 7. Scaling behaviour of the square mean horizon- 
tal width W 2 (I) of the chemical distance for clusters of size 
L — 128 and different values of r\. 




FIG. 4. A schematic representation of an aggregate (left), 
of its backbone (middle) and of its chemical distance (right). 




FIG. 5. The backbone of the cluster shown in Fig. 3. 




FIG. 6. The chemical distance of the cluster of Fig. 3. 
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• • QBM V|=1.0, a=l 


a --□ QBM 11=2.0, a=l 




* * QBM 11 =3.0, a=l 
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FIG. 8. Acceptation profiles for the variables Xi associated 
to the grown bonds, for different values of rj, after t = 500 
time steps. The absence of a threshold value for a(x) implies 
a smooth growth (not by avalanches) . Moreover, one can note 
a dependance on the value of rj, suggesting that in the limit 
r\ = we recover the IP dynamics with a threshold. 



(a) (b) 




FIG. 9. (a) density function of bond 2 at time to; (b) den- 
sity function p2,t +i(3/2) of bond 2 at time to + 1, under the 
hyphotesis of constant stress field. This density concentrates 
more on high values, (c) If one adds the geometric screening 
of the stress field, the P2,t +i(y2) concentrates even more on 
higher values, reflecting the cooperation of the two kinds of 
screening. 
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